NPS55-80-010 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


1 


POSITION  DETERMINATION  WITH 
LORAN-C  TRIPLETS  AND  THE  HEWLETT-PACKARD 
HP-67/9  7  PROGRAMMABLE  CALCULATORS 


by 


R.  H.  Shudde 


March  1980 


Approved  for  public  release;  distribution  unlimited 
Prepared  for: 


f  Naval  Research 
FEDDOCS  =>n'  Virginia  22217 

D  208.14/2:NPS-55-80-010 


KNOXU_BRM« 


MNTCREt  CA  93940 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 


Rear  Admiral  J.  J. 
Superintendent 


Ekelund 


Jack  R.  Borsting 
Provost 


The  work  reported  herein  was  supported  by  the  Office  of  Naval  Research, 
Fleet  Activity  Support  Division,  Code  2  30  and  the  Commander,  Patrol  Wings 
Pacific  as  part  of  the  Tactical  Development  and  Evaluation  Program. 

Reproduction  of  all  or  part  of  this  report  is  authorized. 

This  report  was  prepared  by:  y  rf '  / 


UNCLASSIFIED 


SECURITY   CLASSIFICATION   OF   THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.     REPORT   NUMBER 

NPS55-80-010 


2.  GOVT   ACCESSION  NO 


3.     RECIPIENTS  CATALOG   NUMBER 


4.     TITLE  (and  Subtitle) 

Position  Determination  With  Loran-C  Triplets 
and  the  Hewlett-Packard  HP-67/97  Programmable 
Calculators 


5.     TYPE  OF   REPORT   &   PERIOD  COVERED 

Technical 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.     AUTHORfsJ 

R.    H.    Shudde 


8.  CONTRACT  OR  GRANT  NUMBERfs; 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  CA  93940 


10.     PROGRAM   ELEMENT.  PROJECT,    TASK 
AREA   ft    WORK   UNIT   NUMBERS 

65155N;  R0130,NR  128-217 
N0001479WR90234 


II.     CONTROLLING  OFFICE   NAME   AND   ADDRESS 

Naval  Postgraduate  School 
Monterey,  CA.   93940 


12.     REPORT   DATE 

March    1980 


13.     NUMBER  OF   PAGES 
53 


T4.     MONITORING  AGENCY  NAME  ft    ADDRESSf//  different  from  Controlling  OHIce) 


15.     SECURITY   CLASS,   (of  thia  report) 

Unclassified 


15a.     DECL  ASSIFI  CATION  /DOWN  GRADING 
SCHEDULE 


16.     DISTRIBUTION   ST  AT  EMEN  T  (ol  this  Report) 

Approved  for  public  release;  distribution  unlimited. 


17.     DISTRIBUTION  STATEMENT  (ol  the  abstract  entered  In  Block  20,  It  dllferent  Irom  Report) 


18.     SUPPLEMENTARY   NOTES 


19.     KEY  WORDS  (Continue  on  reverse  aide  It  necessary  and  Identity  by  block  number) 

Loran,  Fixing,  Programmable  Calculator 

Loran-C,  Hyperbolic  Fixing, 

Navigation,  Radio  Positioning, 

Position  Determination,  Calculator 


20.     ABSTRACT  (Continue  on  reverse  aide  It  necessary  and  Identity  by  block  number) 

This  report  presents  an  algorithm  and  HP-67/97  programs  for  position  deter- 
mination with  Loran-C  chains.   Operational  data  cards  are  prepared  in  advance 
for  Loran-C  triplets.   Position  determination  is  performed  using  a  single 
program  card  and  an  appropriate  operational  data  card. 


DD  ,  :°NRM73  1473 


EDITION  OF    1   NOV  65  IS  OBSOLETE 
S/N    0102-014-  6601 


UNCLASSIFIED 


SECURITY  CLASSIFICATION   OF   THIS  PAGE  (When  Data  Entered) 


POSITION  DETERMINATION  WITH  LORAN-C  TRIPLETS 
AND  THE  HEWLETT-PACKARD  HP-6  7/9  7  PROGRAMMABLE  CALCULATORS 


by 


R.  H.  Shudde 

Naval  Postgraduate  School 
Monterey,  California 


March  1980 


The  programs   in  this  report 
are  for  use  within  the  Navy, 
and  they  are  presented  with- 
out representation  or  warranty 
of  any  kind. 


aSJSS 

CONTENTS 

Page 

A.  Introduction   1 

B .  Program  Description   4 

C.  HP-67/97  Calculator  Programs   6 

1 .  User  Instructions   6 

2 .  Sample  Problem   9 

3.  Program  Storage  Allocations  and  Program 

Listings   12 

D.  Loran-C  Fixing  Algorithms 23 

E.  The  Reverse  (Inverse)  Solution  Algorithm   30 

F.  The  Direct  Solution  Algorithm   32 

G.  Discussion  and  Some  Typical  Results   34 

H.     References   38 

APPENDIX:   Loran-C  Station  Parameters   39 


ABSTRACT 

This  report  presents  an  algorithm  and  HP-67/97  programs 
for  position  determination  with  Loran-C  chains.   Operational 
data  cards  are  prepared  in  advance  for  Loran-C  triplets. 
Position  determination  is  performed  using  a  single  program  card 
and  an  appropriate  operational  data  card. 


A.   Introduction 

The  Loran  system  is  a  radio  aid  to  navigation  which 
utilizes  the  principle  of  hyperbolic  fixing.   The  locus  of 
points  for  which  the  difference  in  arrival  time  of  synchronized 
signals  from  a  pair  of  transmitters  is  constant  determines  a 
hyperbolic  line  of  positions  (LOP) .   The  intersection  of  two 
hyperbolic  lines  of  position  from  two  pairs  of  transmitters 
determines  position  or  a  hyperbolic  fix.   That  two  pairs  of 
stations  are  required  for  a  fix  does  not  necessarily  mean  that 
there  are  four  separate  stations,  for  one  station  of  one  pair 
may  be  colocated  with  one  station  of  the  other  pair  forming 
a  Loran    triplet    (Figure  1).   Triplets  may  be  joined  "end-to-end" 
by  station  colocation  to  form  a  Loran    chain     (Figure  2) .   Loran 
chains  are  common  on  both  the  East  and  West  Coasts  of  the 
North  American  continent. 

The  early  "Standard  Loran"  or  Loran-A"  operating  at  a 
frequency  just  below  2MHz  is  still  in  use  in  the  Pacific  area. 
The  present  day  "Loran-C"  operates  at  100-kHz  and  is  in  use 
in  both  the  Atlantic  and  Pacific  Areas.   The  computational 
algorithm  and  programs  described  herein  can  be  used  for  position 
determination  with  Loran-C  triplets.   Further  information  on 
the  history,  development  and  operation  of  the  Loran  systems 
may  be  found  in  References  1  and  2. 


M1'M2 


(a)   Colocated  Master  Stations 


M. 


M, 


(b)   Colocated  Slave  Stations 


M. 


(c)   Colocated  Master  and  Slave 


Figure  1.   Loran  Triplets 
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M  .M 


M  ,M 
u3'  4 


Figure  2.   Loran  Chain  of  Five  Loran  Triplets 


B.   Program  Description 

One  program  card  and  one  operational  data  card 
(described  below)  are  all  that  is  required  for  on-location 
position  determination  from  Loran  triplet  time-difference 
measurements.   Two  program  cards  are  required  to  prepare 
operational  data  cards;  these  operational  data  cards  should 
be  prepared  and  validated  prior  to  on-location  navigational  use. 
Thus  although  three  program  cards  are  described  only  one  program 
card  is  required  for  navigation;  two  program  cards  are  used  to 
prepare  operational  data  cards  during  or  prior  to  mission  planninc 
The  function  of  each  program  card  and  its  intended  use  follows. 


Program  Card  1.  This  program  card  is  used  to  prepare 
master  data  cards.  A  master  data  card  requires  the  following 
information  for  a  master (M)  station/slave (S)  station  pair: 

1.  A  M/S  pair  identification  number. 

2.  The  quantity   At  which  is  the  sum  of  the  coding  delay  plus 
the  one  way  base  line  time  in  microseconds. 

3.  The  latitude  and  longitude  of  the  master  station. 

4.  The  latitude  and  longitude  of  the  slave  station. 

Some  preprocessing  of  these  data  is  performed  before  the  master 
data  card  is  generated.   The  data  generated  require  only  one  side 
of  an  HP-67/97  magnetic  card  for  each  M/S  pair,  thus  a  second 
M/S  pair  may  be  placed  on  side  2  of  the  card  (thus  conserving 
cards)  if  desired.   It  is  envisaged  that  a  master  data  card 
will  be  prepared  in  advance  for  each  M/S  pair  that  might  be 
received  within  an  area  of  operation. 
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Program  Card  2 .   This  program  card  is  used  to  prepare 
an  operational    data    card    for  every  Loran  triplet  within  an 
operational  area.   Each  operational  data  card  contains  data 
merged  from  the  master  data  cards  which  contain  M/S  pair 
information  for  each  pair  of  the  triplet.   These  merged  data 
are  validity  checked,  colocation  of  master  or  slave  determined 
and  encoded. 

The  only  inputs  required  for  this  program  are  the  two 
master  data  cards  that  comprise  the  Loran  triplet.  It  is 
possible  to  prepare  and  store  operational  data  cards  rather 
than  master  data  cards.   This  may  be  desirable  if  there  is  no 
scarcity  of  cards  and  storage  space,  however  the  number  of 
possible  Loran  triplets  is  considerably  larger  than  the 
number  of  M/S  pairs. 

Program  Card  3 .  This  program  card  is  used  in  conjunc- 
tion with  an  operational  data  card  for  position  determination. 
Required  input  is  the  indicated  time  difference   T   for  each 
M/S  pair  of  the  triplet.   Output  is  the  computed  latitude  and 
longitude  of  the  fix.   Note:   Every  Loran  fix  has  two  possible 
solutions.   The  unwanted  solution  can  almost  always  be  rejected 
by  inspection,  however,  if  the  stations  of  the  Loran  triplet 
are  nearly  aligned  then  either  solution  may  be  valid  even 
though  only  one  solution  should  be  consistent  with  the  flight 
plan. 


C.   HP-67/97  Calculator  Programs 
1.   User  Instructions 


CARD  1 


Step 

Instructions 

Input 
Data/Units 

Keys 

Output 
Data/Units 

1. 

2. 

3. 

4a. 
b. 

5a. 
b. 

6a. 

b. 

Read  in  program  card  (both 
sides) 

Input  a  unique  ID  number  for 
the  Loran  pair* 

Input  the  coding  delay   At 

Input  the  master  station  latitude 
and  longitude  (CHS  for  West)** 

Input  the  slave  station  latitude 
and  longitude  (CHS  for  West) 

Run 

Pass  a  blank  data  card  through 

the  card  reader. 

ID 

At 

<J>M 
XM 

xs 

None 

f   a 
f   c 

A 

+ 
C 

E 

ID 

At 

crd 

* 
Note:   Lor an  pairs  are  coded  on  the  navigation  maps  using 

designators  such  as  9930X,  9930Y,  9930Z  and  9930W.   It  is 

suggested  that  the  ID*s  for  these  pairs  be  coded  as 

9930.1,  9930.2,  9930.3  and  9930.4,  respectively.   However, 

any  consistent  scheme  is  acceptable. 

** 

The  format  for  position  data  input  is  of  the  form: 

+DDD.MMSSFF,  where 

DDD  denotes  degrees 
MMM  denotes  minutes 

SS   denotes  seconds 

FF  dnotes  hundredths  of  a  second. 

The  minus  sign  (-)  denotes  Southern  latitudes  or  Western 
longitudes. 

CARD    2 


Step 

Instructions 

Input 
Data/Units 

Keys 

Output 
Data/Units 

1. 

2a. 
b. 

c. 

d. 
e. 

Read  in  program  card  2 

Start 

9.00  will  flash  in  the  display. 
Insert  a  master  data  card 
containing   the  first  pair  of 
a  Loran  triplet. 

9.00  will  flash  in  the  display 
once  more.   Insert  a  master 
data  card  containing  the  second 
pair  of  the  Loran  triplet. 

If  the  data  form  a  proper 
triplet,  "crd"  will  appear  in 
the  display. 

Pass  both  sides  of  a  blank  card 
thru  the  card  reader  to  produce 
the  operational  data  card  for 
the  Loran  triplet. 

E 

9.00 
9.00 

crd 

Should  "error"  appear  in  the  display,  then  the  two  master  data 
cards  do  not  compare  to  form  a  Loran  triplet.   Both  the  latitude 
and  longitude  of  the  colocated  stations  must  be  identical  on 
both  master  data  cards  in  order  to  successfully  produce  an 
operational  data  card. 

Label  the  A  key  position  with  the  identification  number  of  the 
first  Loran  pair  (the  pair  inserted  in  Step  2b)  and  label  the  B 
key  position  with  the  identification  number  of  the  second  Loran 
pair  (from  Step  2c). 

CARD    3 


Step 

Instructions 

Input 
Data/Units 

Keys 

Output 
Data/Units 

1. 

Read  in  both  sides  of  the  program 
card  3. 

2. 

Read  both  sides  of  the  operational 
data  card  for  the  Loran  triplet 
that  you  are  receiving. 

3a. 

Set  to  compute  Solution  A. 

f   a 



b. 

Set  to  compute  Solution  B. 

f  b 



4. 

Input  the  observed  time  delay 
from  the  first  Loran  pair. 

T 

A 



5. 

Input  the  observed  time  delay 
from  the  second  Loran  pair. 

T 

B 



6. 

Compute  fix          Latitude 

Longitude 

E 
R/S 

Latitude 
Longitude 

7. 

Repeat  from  Step  2  with  a  new 
operational  data  card  or  from 
Steps  3  or  4  as  required. 

2 .   Sample  Problem 


CARD  1 


Step 


Instructions 


Input 
Data/Units 


Keys 


Output 
Data/Units 


2. 
3. 

4a. 
b. 

5a. 

b. 

6a. 
b. 

7. 


9a. 

b. 
10a. 

b. 

11a. 
b. 


12. 


In  this  series  of  examples  we 
will  prepare  and  use  data 
cards  for  the  Loran-C  pairs 
9930X  and  9930Y. 

Read  in  program  card  1 
(both  sides) 

Input  the  ID  for  9930X 

Input  the  coding  delay   At 
for  9930X 

Input  the  master  station 

latitude 

and  longitude  (CHS  for  West) 

Input  the  slave  station 

latitude 

and  longitude  (CHS  for  West) 

Compute 

Pass  a  blank  data  card  through 
the  card  reader.   Label  the 
card  9930X  MASTER 

Input  the  ID  for  9930Y 

Input  the  coding  delay   At 
for  9930Y 

Input  the  master  station 

latitude 

and  longitude  (CHS  for  West) 

Input  the  slave  station 

latitude 

and  longitude  (CHS  for  West) 

Compute 

Pass  a  blank  data  card  (or 
the  second  side  of  the  card 
used  in  Step  6b)  through  the 
card  reader.   Label  the  side 
9930Y  MASTER 


These  two  cards  will  be  used 
in  the  next  example. 


9930.1 


36389.66 


34.034604 
■77.544676 


46.463218 
■53.102816 

None 


9930.2 


52541.31 


34.034604 
•77.544676 


41.151193 
■69.583909 

None 


f   a 


f   c 


f   a 


f   c 


f 

A 


9930.10 


36389.66 


crd 


9930.2 


52541.31 


crd 


CARD    2 


Input 

Output 

Step 

Instructions 

Data/Units 

Keys 

Data/Units 

1. 

Read  in  program  card  2. 

2a. 

Start 



E 

9.00 

b. 

While  9.00  is  flashing  in  the 
display,  insert  the  MASTER 
data  card  for  station  9930X 
into  the  card  reader. 

c. 

When  9.00  starts  flashing  in 
the  display  again,  insert 
the  MASTER  data  card  for 
station  9930Y  into  the  card 
reader. 

9.00 

d. 

"Crd'  will  appear  in  the 
display. 



crd 

e . 

Pass  both  sides  of  a  blank 
card  through  the  card 
reader. 

Label  this  card  9930X/9930Y 
OPERATIONAL  DATA  CARD.   Then 
label  the  A  key  position  9990X 
and  the  B  key  position  9930Y. 
This  card  will  be  used  in  the 
next  example. 
. 

0.00 
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CARD    3 


Step 

Instructions 

Input 
Data/Units 

Keys 

Output 
Data/Units 

You  are  receiving  9930X  and 
9930Y  and  wish  to  obtain  a 
fix. 

1. 

Read  in  program  card  3 
(both  sides) 



2. 

Read  in  the  operational  data 
card  for  the  triplet 
9930X/9930Y  (both  sides) 

i 



3. 

Select  Solution  A. 

f   a 

4. 

The  indicated  time  delay  is 

49400  Us  from  9930Y. 

Input  the  indicated  time  delay. 

49400 

B 

0.00 

5. 

The  indicated  time  delay  is 

28800  us  from  9930X. 

Input  the  indicated  time  delay. 

28800 

A 

0.00 

6. 

Solution  A:  42044'57"N  Latitude 
41°07I32"W  Longitude 

f 

Solution  B:  27°00'07"s  Latitude 
102°27,12"E  Longitude 
] Since  you  are  navigating  over    [ 
the  North  Atlantic,  Solution  A 
is  the  desired  fix. 

J 

E 
R/S 

42.4457 
-41.0732 

7. 

Repeat  from  Step  2  with  a  new 
operational  data  card  or  from 
Steps  3  or  4  as  required. 
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Program  Storage  Allocations  and  Program  Listings 
Card  1. 


Registers: 


RO 
Rl 
R2 
R3 
R4 
R5 
R6 
R7 
R8 
R9 


ID 

At 
2c 

M 
M 


MS 


SM 


SO: 

L 

SI: 

T 

S2: 

U 

S3: 

V 

S4: 

X 

S5: 

Y 

S6: 

5    d 

S7: 

AX1 
m 

S8: 

d 

S9: 

f 

RA: 

RB: 

RC: 

e 

m 

RD: 

A6 

] 

RE: 

AA 

RI: 

AX 

m 


m 


Initial  Flag  Status  and  Use: 

0:    OFF,  Unused  2:   OFF,  Unused 

1:    OFF,  Unused  3:   OFF,  Unused 

Display  Status: 

DSP  4,  FIX,  DEC 


User  Control  Keys : 


M     M 


"S    *  xs 


Prepare  data  card 


a  : 
b: 
c: 
d: 
e  : 
12 


Station  ID 


At 


Card  2 


Registers : 


RO 
Rl 
R2 
R3 
R4 
R5 
R6 
R7 
R8 
R9 


+  ID. 


At. 


2c 


A 


Al 
Al 
'Al 
'Bl 
Bl 
'Bl 


SO 
SI 
S2 
S3 
S4 
S5 
S6 
S7 
S8 
S9 


+  ID. 


At. 


2c, 


X 


<; 


A 


A2 
A2 
A2 
B2 
B2 
B2 


RA 
RB 
RC 
RD 
RE 
RI 


a  =  21295.87 
P 

f  =  1/298.26 


Initial  Flag  Status  and  Use: 
0:   OFF,  Vertex  determination 
1:   OFF,  Unused 


2:   OFF,  Validity  checking 
3:   OFF,  Unused 


Display  Status: 

DSP  2,  FIX,  DEG 


User  Control  Keys 

A: 

B 

C 

D 

E:   Run 


a  : 
b: 

c  : 
d: 
e  : 
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Card  3. 

Reg: 

Lsters : 

RO: 

ID1 

Rl: 

At, 

R2: 

(2c)1 

R3: 

A1#  c2,  P 

R4: 

1            F 

R5: 

h  =  V 

R6: 

h 

R7: 

C,,  c±,    H 

R8: 

B1,  S/a  = 

R9: 

a . 

SO: 

ID2 

RA: 

M 

SI: 

At2 

RB: 

u,  N 

S2: 

(2c)2 

RC: 

D,  d 

S3: 

A2 

RD: 

Aa 

S4: 

92 

RE: 

a  =  21295.87 
P 

S5: 

A2 

RF: 

f  =  1/298.26 

S6: 

?2 

S7: 

C2 

S8: 

B2 

S9: 

Initial  Flag  Status  and  Use: 

0:   OFF,  Soln  A,  Soln  B       2:   OFF,  M/S  Vertex  Flag 

1:   OFF,  Unused  3:   OFF,  Unused 

Display  Status: 

DSP  2,  FIX,  DEG 


User  Control  Keys 

A 
B 
C 
D 
E:   Run 


a:   Soln  A 
b:   Soln  B 
c: 
d: 
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c 

0 

■H 

4J 

n 

rH 

O 

w 

H 

rd 

s 

aj 

c 

HI 

Pi 

u 
o 

HJ 

\ 

C 

T3 

■H 

C 

4J 

as 

:• 

•  ■ 

0 

O) 

<D 

to 

-p 

1-1 

3 

0 

c 

04 

4-> 

■H 

b 

w 

rd 

O 

", 

U 

a. 
+ 


CM 


a:: 
I 
n 
CD 


C 

■H 

en 

i 

a 

< 

r\i 

rsi 

\ 

LO 

,-c 

0 

< 

o 

<j 
c 

■H 

to 

+ 

a- 

<■-'] 

M 

C 

■H 

II 


CD                                               CD                              <| 

O 

O 

r-'i  r-_  "-j-  ii  >  cm  T  »r  co  iij  Uj  uj  r-j 

~i   13   O   T   O   CM    -i   'i3   Cj   --1-   — .   •!> 

1                1                                1 
b-j   'O   'O                           b">   lO   'O           lO 

ro  ro  ro                   ro  ro  ro        ro 

*i-  ' u  » ■<  ■*!-  cm  ro  ro  '-i  ro  bj  'u  »-«  T'j  to  *•!■  ■— i  ro  b-j  o  — « 

'    J    -J-   hi   ■—>    -T-   bj   -  <   T   U  j   <-*■   -1-    T    b"j   ro   »-•   -r    hi   U">   O   CM 
II                                                       1                               1                               It 

in  >a  ^o              ^o                    '•o                    '£>                    b"> 
ro  —i  ro              ro                   ro                   ro                   ro 

■—*   l4j   'O 

cm  to  ro 


T5 

u 


-J  -_l  -J  ■+- 


•••J 


-J 


to  r*.  c*; 


a,  co  ■*- •  C'j  ro  m- 
ro  t  "tt-  ■»  T  •» 
<s  cj  '5j  >3  Ci  >a 


O.J 


cj  co  in 

O  -J  -J 

(--  c_.  .  ) 

CO  Cii  u: 


IU  'M 


r  . 
o 
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D.   Loran-C  Fixing  Algorithms 

The  development  of  the  Loran  fixing  algorithms  in  this 
report  is  presented  in  more  detail  in  a  companion  report 
[Ref.  3]  and  will  not  be  repeated  here. 

The  basic  Loran-C  equation  [Ref.  4]  can  be  written  as 


T  =  [Ts  +  p(Ts)]  -  [TM  +  p(TM)]  +  [TB  +  P(TN)]  +  6    (1) 


where 

T   is  the  "indicated  time  difference"  in  microseconds, 
T  ,T_   is  the  distance,  in  microseconds,  from  the  master 
and  the  slave  to  the  receiver,  respectively, 
T    is  the  distance,  in  microseconds,  between  the 

13 

master  and  the  slave, 
6    is  the  assigned  coding  delay,  in  microseconds,  and 
p(T)   is  the  secondary  phase  correction,  in  microseconds, 
for  an  all  sea  water  path  of  length   T. 


The  quantity 


At  =  [TB  +  P(TB) ]  +  6 


is  a  constant  for  each  master/slave  pair.   The  following 
World  Geodetic  System  1972  (WGS  72)  values  have  been  adopted 
for  Loran-C  navigation  [Ref.  4]: 
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v_  =  299792458  meters/second  is  the  velocity  of  light 
in  free  space, 

n   =  1.000338  is  the  index  of  refraction  of  the  surface 
of  the  earth  for  standard  atmosphere  and  100kHz 
electromagentic  waves, 

a   =  6378135.00  meters  is  the  equatorial  radius  of  the 
earth,  and 

f   =  1/298.26  is  the  flattening  factor  (1  -  b/a  ,  where 
b   is  the  polar  radius)  of  the  earth. 

Accurate  formulas  for  computing  the  secondary  phase 
correction   p(T)   are  contained  in  Reference  4,  but  for  use 
with  the  handheld  calculator  the  following  linear  approximation 
[Ref.  3]  will  be  used: 

p(T)  =  a±   +   a2T  , 

where 

a1  =  -0.321, 
and 

a   =   0.000635. 

Using  this  approximation,  it  is  possible  to  solve  Equation  1 

for  the  quantity   T^  -  Tw.   We  find 

S     M 

Ts  -  TM  =  (T  -  At)/(1  +  a2)  .  (2) 


On  the  surface  of  a  sphere  a  hyperbolic  line  of  position 
can  be  represented  by  the  equation  [Ref.  3,  page  175] 


24 


cos  2a  -  cos  2c  ,_. 

tan  r  =  — = ~ ■ — = = (3) 

sin  2c  cos  u)  +  C  sin  2a 


where  the  origin  of  the  coordinate  system  is  at  the  prime  focus 
of  the  spherical  hyperbola,  2c   is  the  spherical  arc  joining 
the  foci,  2a   is  a  constant  for  any  one  LOP,  and   r   and   w 
are  the  spherical  coordinates  of  a  point  on  the  LOP.   If  the 
base  line  of  the  coordinate  system  is  the  arc  joining  the  foci 
then   co   is  the  spherical  polar  angle  from  the  base  line  to  a 
point   P   on  the  LOP  and   r   is  the  spherical  polar  distance 
(or  arc)  from  the  prime  focus  to   P.   Using  the  Loran  system 
we  take    c,    =  +1   if  the  prime  focus  is  at  a  master  station 
and  we  take  C  =  -1   if  the  prime  focus  is  at  a  slave  station. 

If  we  take   v  =  v„/n   to  be  the  velocity  of  100kHz  electro- 
magnetic radiation  of  the  earth's  surface  then 


or,  using  Eq .  (2) , 


where 


2a  =  v(Ts  -  TM)/ae  , 


2a  =  (T  -  At) /a   ,  (4) 

p 


a  (1  +  a  ) 

a   =  — t —  =  21295.87  ys  . 

P       vQ/n 

The  baseline  between  master  and  slave  can  be  obtained  from 


2c  =  v  TD/a   .  (5) 

B   e 

Here   2c   is  computed  by  program  card  1  (preparation  of  master 
data  cards)  using  the  algorithm  in  Section  E. 
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Consider  a  Loran-C  triplet  with  master  stations  colocated. 
Let   C-i   and   5?   denote  the  azimuth  angles  of  slave  1  (S,  ) 
and  slave  2  (S») ,  respectively,  measured  from  North  toward  the 
East  from  the  master  stations  (M)  (see  Fig.  3) .   Further,  let 
a      and   r   denote  the  azimuth  and  spherical  polar  arc  (distance! 
of  the  receiver  (R)  from  M.   For  this  geometry,  Eq.  (3)  can 
be  written  as 

B. 

tan  r.  =  ~ ; = — r — ; — - —  (6) 

i    C.  cos (a-C • )  +  A. 


where 


and 


A.  =  £ .  sin  2a . 
li        l 


B.  =  cos  2a.  -  cos  2c. 
ill 


C .  =  sin  2c . 
i  l 


for  the  i —  Loran  pair,  i  =  1,2.  Since  r  =  r,  =  r?,  tan  r. 
can  be  eliminated  in  Eq.  (6) .  The  resulting  equation  can  be 
rewritten  as 

C  cos  a  +  S  sin  a  =  K  ,  (7) 


where 


C  =  B,C2  cos  C2  ~  B2Cl  cos  ^l' 
S  =  B,C2  sin  C2  "  B2Cl  sin  ^l' 


and      k  =  BpA,  -  B..A,, 
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R 


Figure  3.   Geometry  of  a  Loran  Triplet  and  a  Receiver 
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If  we  define   p  >  0   and  y     by  the  equations 

p  cos  y  =  C, 
and  (8) 

p  sin  y    =   S, 

then 

p  =  /c2  +  S2  , 

and 

Y  =  qatn(S,C)  . 

Here  the  function   qatn(y,x)   is  the  arctangent  of   y/x 
adjusted  for  the  proper  quadrant  according  to  the  signs  of   x 
and   y.   A  compact  form  of  this  function  is 


qatn(y,x)  =  tan    -^ +  iTt(x  <  0?) 

x  +  10   t(x  =  0?) 


where 

t(z)  =  1   when   z   is  true 
and 

t ( z)  =0   when   z   is  false. 

When  convenient  we  will  use  the  notation   qatn(y/x)   inter- 
changeably with  qatn(y,x) .   The   qatn   function  is  equivalent 
to  the  polar  angle  obtained  using  the  rectangular  to  polar 
conversion  function  on  the  HP-67/97. 

Now  substitute  Eq .  (8)  into  Eq .  (7)  and  solve  for 

a  =  y  +  cos   (k/p)  (9) 
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to  obtain  the  azimuth  angle   a   of  the  two  points  of  intersection 
of  the  LOP's.   Finally  we  obtain  a  value  for   r   by  substituting 
each   a  into  Eq.  (5) .   We  find  that 


r  =  qatn 


B. 

1 


C.  cos (a-C . )  +  A. 
1         l      l 


for   i  =  1  or  2 


The  distance  and  azimuth  from  M  or  the  triplet  vertex  can  be 
converted  into  the  latitude  and  longitude  of  the  two  possible 
positions  of   R. 

The  fixing  algorithm  then  uses   a   and   r   in  the  direct 
solution  algorithm  of  spheroidal  geodesy  (Section  F) . 
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E.   The  Reverse  (Inverse)  Solution  Algorithm 

This  reverse    solution  algorithm  is  a  modification  of  the 

first  order  in  flattening  (f)  algorithm  given  by  Thomas  [Ref.  5, 

pp.  8-10].   Thomas'  notation  has  been  followed  as  closely  as 

possible  for  ease  of  comparison  of  the  algorithms.   The   qatn 

function  is  defined  in  Section  D.   West  longitudes  (A)  and 

South  latitudes  (<p)    are  negative.   We  are  given  the  points 

P  ((j),  ,  A,  )  ,  P„ ($„ , A~)   on  the  spheroid  and  are  to  find  the 

distance   S   between  the  points  and  the  forward  and  back 

azimuths,  a,  ~   and   a?,  .   Given  quantities  are   (f),  ,  A,,  <j>_ 

and   A„.   No  assumptions  about  the  relative  location  of   P, 

and   P9   are  required.   The  modified  reverse    solution  algorithm 

is  : 

6.    =   tan"    [(1-f)    tan   &,.],         i   =    1,    2, 
l  l 

e    =  (e    +ej/2,    A6    =  (e  -en)/2,     aa  =  a0  -  a,, 

ml2  m21  21 

2  2  2  2 

AA   =  AA/2,   H  =  cos  A0   -  sin   9   =  cos   9   -  sin   A9   =  cos  9   cos  0  . 
m  m         m         m  m         1       2 

2  2  2  2 

L  =  sin  A9   +  H  sin  AA   =  sin  (d/2) ,  1  -  L  =  cos  (d/2) , 
m  m 

-1  2        2 

d  =  cos   (1  -  2L) ,  U  =  2  sin   9   cos   £9  /(l  -  L), 

m        m 

2         2 
V  =  2  sin   A9   cos   9  /L,  X=U+V,  Y=U-V, 
m       m 

T  =  d/sin  d,   6nd  =  f(TX  -  Y)/4,  S  =  a  (T-6,d)  sin  d, 

1  e     1 

F  =  2[Y  -  (1  -  2L)  (4  -  X)  ]  ,    G  =  f  T/2 , 

Q  =  -(FG  tan  AA)/4,   AA  '  =  (AA  +  Q) /2 

m 
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t   =  qatn(-sin  A6   cos  AA  '  ,  cos  9   sin  AA'), 
1  m       m       m       m 


t„  =  qatn(cos  A6   cos  AA',  sin  9   sin  AA ' )  , 
2  m       m       m       m 


a!2  =  tl    +    V   a21  =  tl    ~    t2 


This  reverse  solution  algorithm  is  used  by  program  card  1 

(preparation  of  master  data  cards)  to  compute  the  baseline 

distance   2c   and  the  azimuths   C„o   and   L„  between  the 

MS         SM 

master  and  slave  stations  of  a  Loran  pair. 

Details  of  the  modifications  made  to  Thomas'  algorithm 
are  contained  in  Reference  3. 
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F.   The  Direct  Solution  Algorithm 

This  direct    solution  algorithm  is  a  modification  of 
the  first  order  in  flattening  (f)  algorithm  given  by  Thomas 
[Ref.  5,  pp.  7-8] .   Thomas'  notation  has  been  followed  as 
closely  as  possible  for  ease  of  comparison  of  the  algorithms. 
The   qatn   function  is  defined  in  Section  D.   West  longitudes 
and  South  latitudes  are  negative.   We  are  given  the  point 
P  ((f>  ,  A  )   on  the  spherioid,  where   <!>  ,  ^,   are  the  geodetic 
latitude  and  longtiude  (geographic  coordinates) ;  the  forward 
azimuth   a,2   and  the  distance   S   to  a  second  point   P_(4>_,A  ); 
and  from  these  we  are  to  find  the  geographic  coordinates 
^  ,  \        and  the  back  azimuth   a2i'   Tne  given  quantities  are 
$  ,  A  ,  a  „   and   S.   No  assumptions  about  the  relative  loca- 

_1_       X.  -L  « 

tion  of   P,   and   P~   are  required.   The  modified  direct 
solution  algorithm  is: 

6  =  tan   [(1-f)  tan  <J)  ]  ,    M  =  cos  6  sin  a 

2 
N  =  cos  6   cos  a      c   =  fM,   c   =  f(l  -  M  )/4, 

D  =  1  -  2c   -  cM,    P  =  c  /D,    0      =   qatn(N,  sin  9  ) 

d  =  S/(a  D) ,    u  =  2 (Q  -d) ,    W  =  1  -  2P  cos  u, 
e  1 

V  =  cos(u  +  d) ,    Y  =  2PVW  sin  d,    Aa  =  d-Y , 

a   =  qatn[-M,  -  (N  cos  Aa  -  sin  9   sin  Aa ) 1  , 

2  1/2 
K  =  (1-f)  [M  +  (N  cos  Aa  -  sin  6   sin  Aa )  ] 

4>   =  tan"  [(sin  6   cos  Aa  +  N  sin  Aa)/K]  , 

=  qatn (sin  Aa  sin  a.  cos  9n  cos  Aa  -  sin  9   sin  Aa  cos  a   ) 

iz  1  1  ^2 

H  =  c1Aa,   AA  =  An  -  H,    A2  =  \±  +   AA. 
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This  direct  solution  algorithm  is  used  by  program 
card  3  (improved  fix  program)  to  compute  the  latitude  and 
longitude  of  the  receiver  using  the  azimuth  and  range  of  the 
receiver  from  the  Loran  triplet  vertex. 

Details  of  the  modifications  made  to  Thomas'  algorithm 
are  contained  in  Reference  3 . 
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G.       Discussion   and   Some   Typical    Results 

The   HP-67    program  design    specifications    of 
COMPATWINGSPAC    [Ref.    6]    are   contained    in   the    following 
statement. 

"There    is   a  need   for   an   HP-67   program  that   will   compute   a 
geographical  position    from  two   Loran   delay   rate   readings. 
Several   methodologies   are   available   to   compute   the 
desired   position  but   computational   complexities    increase 
with  the   desired  accuracy  and   flexibility.      The  most 
desirable   accuracy   would  be   an   error   of    less   than   4   n.mi. 
at   a   range   of   500   n.mi.    with    less   error   closer   to  the 
stations.       It    is    likely   that  program   length   considerations 
will   require   that   the    station  pairs   have    a   common   site 
(i.e.    two   slaves   or   two  masters   at   the   same    location). 
This    is   not   an   unusual    situation   as   evidenced  by   strings 
of   station  pairs   along   coast    lines.      A  data  card  will 
probably  be   necessary    for   the   station  pairs    to  be   used. 
However,    more   than   one   program  card   is   unacceptable   due 
to   the   decrease    in    functional    utility  when   compared   to 
the   manual   plotting  method.      As   a    final    requirement,    the 
fix   should  be   obtainable   on   either   side   of   the   baselines 
connecting   the    stations,    and   not    limited   to   a   geometric 
position    relative   to   one   side   or    the   other   of   the 
stations.  " 

It  was    further    stated    that    the   maximum  computation    time    to 
obtain    a    fix   be    1.5    minutes. 

It    is    felt    that    these   design   goals   have   been    satisfied 
Although   one    program    is    required    to   prepare   master    data   cards 
for   all    Loran-C   pairs    and   a    second   card    is    required    to   prepare 
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operational  data  cards,  one  each  for  every  triplet,  this  pre- 
paration should  be  done  only  once.   The  data  cards  should  be 
supplied  to  users  verified  and  labeled,  by  the  Fleet  Mission 
Program  Library.   One  program  card  and  an  appropriate  oper- 
ational data  card  are  all  that  is  required  for  the  fixing 
algorithm. 

The  fixing  algorithm  will  display  one  of  the  two 
possible  receiver  positions  in  38  seconds  following  the  entry 
of  the  time  delay  readings.   Since  there  are  situations  in  which 
either    of  the  two  solutions  could  be  the  valid  solution;  the 
decision  of  which  solution  to  use  should  be  left  to  the  operator, 
not  the  program  designer. 

Testing  of  the  algorithm  for  all  Loran-C  triplets  and 
positions  relative  to  those  triplets  was  too  extensive  a 
program  to  be  carried  out  in  the  available  time.   Some  "typical" 
scenarios  however  are  presented  in  Tables  I  through  IV.   As  can 
be  seen  all  errors  are  all  well  within  the  design  specifications 
of  4  n.mi  at   500  n.mi  range  from  the  stations.   The  time  delay 
values  in  these  Tables  were  generated  using  a  program  discussed 
in  Reference  3.   It  is  recommended  that  the  P-3  community  test 
the  algorithm  for  accuracy  in  known  areas  of  operation  and 
examine  the  results  for  possible  regions  in  which  the  algorithm 
may  fall  outside  the  design  requirements.   Such  testing  should 
be  compatible  with  the  known  "unreliable  regions"  shown  on  the 
Loran-C  charts. 
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Table  I.   Moffett  Field  South 


Position 

Indicated 

Time  Delay 

Fix 

Error 

Lat 

Long 

9940X 

9940Y 

Lat(N) 

Long(W) 

n.mi 

24°N 

122°W 

27726.19 

40912.76 

23°59'55" 

122°00'01" 

0.08 

26 

122 

27715.97 

40998.39 

25°59'57" 

122°00,01" 

0.05 

28 

122 

27702.41 

41117.84 

27°59,59" 

122°00'00" 

0.02 

30 

122 

27683.53 

41291.85 

29059.59" 

122°00'00" 

0.02 

32 

122 

27655.47 

41555.46 

32°00,00" 

122"00,00" 

0.00 

34 

122 

27609.63 

41959.57 

34°00'00" 

122°00'00" 

0.00 

36 

122 

27523.56 

42544.11 

36°00'00" 

121059'59" 

0.01 

38 

122 

27334.61 

43248.22 

38°00'00" 

121059'58" 

0.03 

Table  II 


Moffett  Field  West 


Position 

Indicated 

Time  Delay 

Fix 

Error 

Lat 

Long 

9940Y 

9940W 

Lat(N) 

Long(W) 

n.mi 

37°N 

122°W 

42892.86 

16257.23 

36°59'59" 

122°00,01" 

0.02 

37 

125 

43056.68 

15765.13 

37°00,00n 

125o00'00" 

0.00 

37 

128 

43137.78 

15327.12 

37°00'00" 

128o00'00" 

0.00 

37 

131 

43191.10 

14970.77 

37°00'00" 

131°00'00" 

0.00 

37 

134 

43232.38 

14683.74 

37°00'00" 

134°00,00M 

0.00 

37 

137 

43267.42 

14449.40 

37o00,00" 

137°00'00" 

0.00 

37 

140 

43298.80 

14254.02 

37°00'00" 

140°00'01" 

0.01 

37 

143 

43327.85 

14087.43 

37°00'01" 

142°59'59" 

0.02 
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Table  III.   Brunswick  Northeast 


Posi 

tion 

Indicated 

Time  Delay 

Fix 

Error 

Lat 

Long 

7930Z 

9930X 

Lat(N) 

Long (W) 

n.mi 

60°N 

30°W 

52437.86 

28451. 72 

60°00'03" 

29059- 32" 

0.24 

58 

35 

51960.93 

28391. 50 

58°00'00" 

34059. 46" 

0.11 

56 

40 

50992.37 

28359. 15 

55°59'59" 

39059.54.. 

0.06 

54 

45 

49292.46 

28370.85 

53059'59" 

44°59'57" 

0.03 

52 

50 

47165.60 

28490. 64 

52o00'00" 

49°59'59" 

0.01 

50 

55 

45236.59 

29070. 48 

50°00'00" 

55°00'00" 

0.00 

48 

60 

44505.60 

30991. 94 

48o00'00" 

60°00'00" 

0.00 

46 

65 

44475.70 

33697. 14 

46°00'00" 

65°00'00" 

0.00 

44 

70 

44588.91 

36567. 42 

43059.59- 

69°59'59" 

0.02 

Table  IV.   Jacksonville  Southeast 


Position 

Indicated 

Time  Delay 

Fix 

Error 

Lat 

Long 

9930W 

9930X 

Lat (N) 

Long(W) 

n.mi 

9°N 

47°W 

13058.04 

36466.  46 

8°59'19" 

46°59'22" 

0.92 

12 

52 

12984.71 

37288.  35 

11°59'34" 

51°59'37" 

0.57 

15 

57 

12898.73 

38267. 58 

14°59,44" 

56059'47" 

0.34 

18 

62 

12793.91 

39431.  32 

17059'52" 

61°59,54" 

0.16 

21 

67 

12656.52 

40794.  36 

20°59'56" 

66°59'57" 

0.08 

24 

72 

12451.30 

42330.  55 

23°59'59" 

71059'59" 

0.02 

27 

77 

12097.12 

43876. 62 

27o00'01" 

77°00'00" 

0.02 

30 

82 

12973.95 

44768.  53 

30°00'01" 

82°00,06" 

0.09 
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APPENDIX.   Loran-C  Station  Parameters 


The  following  list  contains  the  pertinent  parameters 
for  each  Loran-C  station  pair.   This  list  was  compiled  from 
data  in  Reference  4.   Each  column  contains  the  following 
information: 

1.  The  Loran-C  station  pair  designator 

2.  At,  the  sum  of  the  coding  delay  plus  one  way  baseline 
time,  including  the  secondary  phase  correction  for  an 
all  seawater  path,  in  microseconds. 

3.  The  master  station  latitude. 

4.  The  master  station  longitude. 

5.  The  slave  station  latitude. 

6.  The  slave  station  longitude. 

In  this  list,  negative  longitudes  are  West  and  positive 

longitudes  are  East.   If  desired,  this  convention  may  be 

reversed  since  the  algorithms  are  independent  of  such  external 

conventions;  if  this  is  done,  care  should  be  taken  that  the 

signs  of  all  longitudes  in  the  list  are  reversed.   In  columns 

3  through  6  the  latitudes  and  longitudes  appear  to  be  in 

decimal  form,  but  the  actual  format  is  DDD.MMSSFF  (which  is 

compatible  with  the  HP-67/97  H.MS  input  mode)  where 

DDD  designates  degrees, 

MM  designates  minutes, 

SS  designates  seconds,  and 

FF  designates  hundredths  of  seconds. 
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4390K..  15972.  23 

. 16.444335. 

-169.303120. 

20 

.144916. 

-155.530978 

4399V, 34253. IS. 

16.444395. 

-163.303120. 

IsLC' 

,234177. 

-178. 173820 

593UX, 13131.38 

46.482720. 

-067.553771. 

41 

.  151133. 

-069. 5:3:3909 

5330V, 23755. 02. 

46. 432720. 

-067.553771. 

46. 

463218. 

-053. 182316 

5390X , 1 3343 . 60 

cr  -i   cr -r  cr .-. -? .-. 
■J  1  .  _T  JOl  O. 

-122.220224. 

55 

. 262085 

.-131.151365 

5990V , 23927 . 36 . 

cr  1   cr  -?  cr .-.  ~f .-. 
•J  1  .  J  1  -JO  1  O 

-122.220224. 

47 

. 834799. 

.-119.443953 

5390Z, 42266. 63. 

cr  h   cr  — ■  cr .-.  -?  .-. 
_'l  .  .J  r  JO  r  O. 

-122.220224. 

50 

. 362972 

-127.212935 

7330N, 15063.02 

53.591727. 

-045. 102747. 

64 

. 542653 

.-023.552175 

7930.y...  27303.  77. 

59.591727. 

-045. 102747. 

62 

. 175968. 

-007. 842671 

73302,43212.20. 

59.591727. 

-045. 102747. 

46 

.463213 

-053. 182816 

796 9X, 13304. 45. 

63. 194231. 

-142.433190 

57 

.262021 

.-152.221122 

7960V, 29651. 14. 

63. 194231. 

-142.433190. 

55 

. 262035 

.-131. 151365 

, ■■'370N,  30965.  64. 

62. 175363. 

-007.042671 

cr  A 
■_'-t 

. 432980 

.+803. 173633 

7370X, 15043. 19. 

62. 175363. 

-007.342671 

68 

. 3886 1 5 

.+814.274700 

7970V, 48944. 53. 

62. 175963. 

-007.042671 

•  64 

. 542653 

.-023.552175 

7970Z.I  63216.  38 

62. 175363 

-007. 342671 

70 

.545261 

. -088. 435863 

7988W, 12389.54. 

30.533374. 

-085. 100330 

30 

. 433382 

. -830. 494360 

79S0>:..  27443.38. 

30.533374. 

-yy5.  1 00 930. 

2  b 

.315501 

. -037. 500009 

iV, -45281.  88. 

33.533374. 

-085. 100938 

C-  i 

.015349 

. -888. 065352 

"7  •  j  o  ~i "?  u-r  -i  ="  -i  -  ~»  "' 

•-;  L-i  cr  q  ■-.  -.  -7  4 

-035. 100330. 

34 

. 034604 

.-077.544676 

.  JO'Jii.  .'  •-  i  -!ktc  .  ,  j..  . 

•_■  0  .  J-  y  ■-■  Cm  4  . 

7330;::,',  12755.  97. 

33.522061 

016. 430596. 

■-.cr 

.312033 

.    812.312396 

7990V, 32273. 30. 

33.522361. 

015. 430536. 

40 

. 532835 

.    027.520152 

7  ?  '::|  0  7 ,  5  0  9  9  9 .  69. 

33.522061. 

016. 480596 

42 

. 033643 

.    083. 121590 

14355! 11. 

39.519754. 

-087.291214 

33 

. 593374 

. -085. 188938 

3!y!,  31 162.  06. 

33.510754. 

-337.291214. 

42 

. 425860 

.-076.433386 

3970V. 47753. 74. 

39.510754. 

-037.231214 

48 

. 364334 

.-834.331847 

13695.51. 

34. 0346Q4. 

-077.544676. 

d  1 

.815343 

. -080. 065352 

389. 66. 

34.334604 

-077. 544676 

46 

.463213 

. -053. 102816 

"  ' , '    c;  "■  ="  A  -i     ■"'  i 

34. 334604. 

-077.544676 

41 

. 151133 

-863. 583989 

"??.   £3560.72. 

34. 334604. 

-077. 544676. 

39 

.518754 

.-837.291214 

~i\, j ,  1 3736 .  99 

33. 330662 

-113.495637. 

47 

. 034739 

■-119.443953 

":i4i*'-  23094,50 

39. 330662 

-113.435637. 

•~i  i*j 

, 465639 

.-122.234453 

jV . 41957, 33 

33. 330662. 

-113.435637. 

35 

. 191813 

.-114.481743 

9363N. 13797.20. 

42. 425060. 

-076.433336. 

46 

. 432728 

-867.553771 

-  '"'  '•"  "1   ,  f:  "?6  9 .  94; 

42. 425060. 

-F17S,  4  Q :";:"-: >?,fk 

41 

, 151133. 

-069. 583909 

3963  V . 42221 . 65. 

42. 425060. 

-076. 493336. 

34 

, 094604 

-877.544676 

0Z, 57162. 06. 

42.425063. 

-076. 433386. 

39 

510754. 

-037.291214 

3970H, 15233. 94. 

24.43041    . 

141. 13290    . 

24 

17077    . 

153.58515 

j  :i  ~>  i~i  |  /  ~  ;~  |T  c,  k:  |  ~, 

24.43041    . 

141. 19298    . 

42 

443700. 

143. 438986 

~!970V  .  59463.  1  3. 

24.43041    . 

141. 19290    . 

26 

, 36^439. 

128.085621 

99737, 30746. 73. 

24.43341    . 

141.1 9290    . 

09 

.324566. 

138.035523 

3330K  ■  l-i  375.  32. 

57. 090933. 

-170.145931. 

52 

. 434505 

.+173.185231 

3393V ■ 32069 « 03. 

57.  0909:39 

-170. 145931. 

65 

.  144012. 

-166.531447 

3380Z  •  46593.  19. 

57. 090938. 

-170. 145931. 

cr-? 
•-.'  1 

, 26202 1 

-152.221122 

40 


Coverage  of  Loran-C  Systems 

Station  Location 

4990  Central  Pacific 

5930  East  Coast,  Canada 

5990  West  Coast,  Canada 

7930  North  Atlantic 

7960  Gulf  of  Alaska 

7970  Norwegian  Sea 

7980  Southeast  U.S.A. 

7990  Mediterranean  Sea 

8970  Great  Lakes 

9930  East  Coast,  U.S.A. 

9940  West  Coast,  U.S.A. 

9960  Northeast  U.S.A. 

9970  Northwest  Pacific 

9990  North  Pacific 
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INITIAL  DISTRIBUTION  LIST 


NO.  OF  COPIES 


Defense  Technical  Information  Center  12 

Cameron  Station 
Alexandria,  VA   22314 

Library,  Code  0212  2 

Naval  Postgraduate  School 
Monterey,  CA   9  39  40 

Library,  Code  55  2 

Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA   93940 

Dean  of  Research,  Code  012  1 

Naval  Postgraduate  School 
Monterey,  CA   9  3940 

Office  of  Naval  Research  2 

Fleet  Activity  Support  Division 

Code  ONR-2  30 

800  North  Quincy  Street 

Arlington,  VA   22217 

Attn:   Mr.  Robert  Miller 

Office  of  Naval  Research  1 

Code  ONR-4  34 

80  0  North  Quincy  Street 

Arlington,  VA   22217 

Navy  Tactical  Support  Activity  2 

P.O.  Box  1042 

Silver  SDrings,  MD   20910 

COMPATW1NGSPAC  4 

Naval  Air  Station 
Moffett  Field,  CA   94035 
Attn:  Code  51  and  Code  5  32 

COMPATW1.MGSLANT  2 

Naval  Air  Station 
Brunswick,  ME   04011 
Attn:   Code  N7 

COMPATWING  ELEVEN  TSC  NAS  1 

Jacksonville,  PL   322.12 

COMSEABASEDASWW1NCSLANT  2 

Nava]  Air  Station 
Jackson v  i I  lo ,  FL   32223 
Attn:  TDi.E  Officer 
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NO.  OF  COPIES 


Commanding  Officer  1 

Air  Test  and  Evaluation  Squadron  1  (VX-1) 
Patuxent  River,  MD   20670 
Attn:  Code  713 

Commanding  Officer  1 

Submarine  Development  Group  Two 
Groton,  CT   06340 

Director  1 

Strategic  Systems  Project  Office 
1931  Jefferson  Davis  Highway 
Arlington,  VA   20376 
Attn:  Code  SP2021 

Naval  Air  Development  Center  1 

Johnsville,  PA   18974 
Attn:  Code  2022 

Center  for  Naval  Analysis  1 

1401  Wilson  Boulevard 
Arlington,  VA   22209 
Attn:   Greg  Watson 

Naval  Weapons  Laboratory  1 

Dahlgren,  VA   22443 

Naval  Weapons  Center  1 

China  Lake,  CA   93555 

Naval  Surface  Weapons  Center  1 

White  Oak 

Silver   Spring,  MD   20910 

Naval  Research  Laboratory  1 

Washington,  D.C.   20390 

David  Taylor  Naval  Ship  Research  1 

and  Development  Center 
Bethesda,  MD   20034 

Naval  Ocean  Systems  Center  1 

San  Diego,  CA   92132 

Naval  Intelligence  Support  Center  1 

4301  Suitland  Road 
Washington,  D.C.   20390 

Naval  Electronics  Systems  Command  1 

2  511  Jefferson  Davis  Highway 
Arlington,  VA   20  360 
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NO.  OF  COPIES 


Naval  Underwater  Systems  Center  1 

Code  SA3  3 

New  London,  CT   06320 

Naval  Ship  Engineering  Center  1 

Hyattsville,  MD   20782 

Naval  Coastal  Systems  Laboratory  1 

Panama  City,  FL   32401 

Naval  Air  Systems  Command  1 

Code  3  70 

Washington,  D.C.   20361 

Naval  Sea  Systems  Command  1 

Code  0  3424 

Washington,  D.C.   20362 

Naval  Underwater  Systems  Center  1 

Newport,  Rhode  Isla:id   02840 

Naval  Ordnance  Station  1 

Indian  Head,  MD   2  0 S40 

Naval  Surface  Weapo  is  Center  1 

Dahlgren,  VA   224  48 

Anti-Submarine  Warfare  Systems  Project  Office  1 

Department  of  the  Navy 
Washington,  D.C.   20360 
Attn:   Code  ASW-137 

Surface  Warfare  Development  Group  2 

Naval  Amphibious  Base 

Little  Creek,  VA   23511 

Attn:   Code  N32  and  C.  Cartledge 

COMNAVSURFLANT  1 

Norfolk,  FA   23511 
Attn:  Code  N612 

Deputy  Commander  1 

Operational  Test  and  Evaluation  Force  Pacific 
NAS  North  Island 
San  Diego,  CA   92  135 
Attn:  Code  202 

VP-24  1 

NAS 

Jacksonville,  FL   32212 

Attn:   LTJG  Kenneth  D.  Walker 
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Assistant  Wing  TAC  D&E  Officer  1 

ASW  Operations  Department 

NAS 

Cecil  Field,  FL   32215 

Attn:   Code  VS-24 

Naval  Electromagnetic  Spectrum  Center  1 

Naval  Communications  Unit 
Washington,  D.C.   20390 
Attn:   CDR  Claude  JaVarre 

VS-Support  Unit  1 

NAS 

Cecil  Field,  FL   32215 

Attn:   CWO  Wentworth 

PROCAL-Code  18  1 

Fleet  Readiness  Office 

NOSC 

San  Diego,  CA   92152 

Attn:   Jim  Grant 

Chief  of  Naval  Materiel  1 

MAT  08D1 

Department  of  the  Navy 

Washington,  D.C.   20360 

Attn:   CDR  Philip  Harvey 

LT  Jan  Smith  1 

COMCRUDESGRU  TWO 

FPO  New  York,  NY   09501 

Dr.  Thomas  Burnett  1 

Dept .  of  Management  Science 
OR  Study  Group 
U.S.  Naval  Academy 
Annapolis,  MD   21401 

U.  S.  Naval  Oceanographic  Office  1 

Washington,  D.C.   20309 
Attn:   Paul  D.  Thomas 

Navigation  Department  I 

Defense  Mapping  Agency 
Hydrographic/Topographic  Center 
Washington,  D.C.   20315 
Attn:   G.  R.  Deyoung 
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LT  Michael  D.  Redshaw  1 

Helicopter  Antisubmarine  Squadron  7 
FPO  New  York   09501 

Mr.  Raymond  F.  Fish  1 

Naval  Underwater  Systems  Center 
Newport,  RI   02840 

LT  Kenneth  W.  Peters  1 

Commanding  Officer 

U.  S.  Naval  Air  Facility 

APO  New  York   09  40  6 

LT  Peter  W.  Marzluff  1 

Surface  Warfare  Officer's  School 
Newport,  RI   02840 

LT  Ray  C.  Pilcher,  Jr.  1 

Surface  Warfare  Officer's  School 
Newport,  RI   0  28  40 

LT  William  M.  Yerkes  1 

ASWOC  NAF 
S.^gonella  Sicily 
FPO  New  York  09523 

LT  A.  J.  Koorey  1 

Commander 

Cruiser-Destroyer  Group  Two 
FPO  New  York   09501 

LCDR  E.  G.  Schwier  1 

USS  Estocin  (FFG  15) 
FPO  New  York   09501 

LCDR  R.  R.  Hillyer  1 

Code  N311 

Commander,  Carrier  Group  Four 

FPO  New  York   09501 

LCDR  J.  W.  Pattison  1 

Code  N321 

Commander,  Carrier  Group  Eight 

FPO  New  York   09501 

LT  Richard  C.  Myers  1 

ASWOC  NAF  LAJES 

Azores 

APO  New  York   09523 
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LT  R.  N.  Christianson  1 

ASW 

Naval  Postgraduate  School 

Monterey,  CA   93940 

LT  L.  R.  Erazo  1 

ASWOC  SIGONELLA 
FPO  New  York   09523 

LT  R.  E.  Springman  1 

USS  Lexington  (ANT-16) 
Naval  Air  Station 
Pensacola,  FL  32504 

LT  Philip  D.  Ward  1 

Patrol  Squadron  26 
Special  Project  Det 
Naval  Air  Station 
Brunswick,  ME   04011 

LT  Robert  J.  Knight  1 

USS  John  F.  Kennedy  (CV-67) 
FPO  New  York   09501 

LT  Richard  0.  White,  Jr.  1 

USS  Dwight  D.  Eisenhower  (CVN-69) 
FPO  New  York   09501 

LT  David  E.  Hebdon  1 

USS  John  F.  Kennedy  (CU-67) 
FPO  New  York   09501 

LCDR  Michael  A.  Santoro  1 

Air  Antisubron  41 
NAS  North  Island 
San  Diego,  CA   92135 

LCDR  Robert  M.  Hanson  1 

Helasron  Four 
NAS  North  Island 
San  Diego,  CA   92135 

LT  Michael  D.  Thomas  1 

2476  Ridgecrest  Ave. 
Orange  Park,  FL   32073 

LT  Carl  E.  Garrett,  Jr.  1 

Class  65,  Dept.  Head  Course 
Surface  Warfare  School,  Naval  Base 
Newport,  RI   0  2840 
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Naval  Postgraduate  School 
Monterey,  CA   93940 

Attn:   Prof.  R.  N.  Forrest,  Code  55Fo  1 

Prof.  R.  H.  Shudde,  Code  55Su  200 

R.  J.  Stampfel,  Code  55  1 


LT  M.  D.  Clary 

USS  PAIUTE  (ATF-159) 

FPO  New  York  09501 
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